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We introduce a real-space exact renormalization group method to find exactly solvable quantum 
spin chains and their ground states. This method allows us to provide a complete list for exact 
solutions within SU(2) symmetric quantum spin chains with S < 4 and nearest-neighbor interactions, 
as well as examples with 5* = 5. We obtain two classes of solutions: One of them converges to the 
fixed points of renormalization group and the ground states are matrix product states. Another one 
does not have renormalization fixed points and the ground states are partially ferromagnetic states. 
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I. INTRODUCTION 

Understanding the physical properties of quantum 
many-body systems is an important common issue in 
condensed matter physics and quantum information the- 
ory. The number of parameters required to describe a 
random state grows exponentially with the number of 
particles, which makes the computation of many-body 
systems very difficult, even numerically [1 . However, re- 
cent development in quantum information theory implies 
that only a corner of such a huge Hilbert space is relevant 
for describing the low-energy states of physical systems 
[2j [3] . The characteristic feature of this corner seems to 
be an area law [4 : the von Neumann entropy of a subsys- 
tem in the many-body ground state scales with the bor- 
der area, rather than the volume - the case for a random 
state. This means that the ground states of quantum 
many-body systems usually only contain a small amount 
of entanglement. It is natural to take this advantage and 
design clever parametrizations of states which both cap- 
ture the essential physics and allow classical simulations 
with a polynomial time. In one dimension (ID), the ma- 
trix product state (MPS) [5]-[7] is a candidate for such 
a purpose. The MPS lies at the heart of the success of 
the density-matrix renormalization group (DMRG) [8]- 
[12], which has been proved to be an accurate numerical 
method for describing the low-energy states of quantum 
lattice models. Recently, there are many interesting ex- 
tensions along this direction, including infinite MPS for 
critical systems [13 , continuous MPS for quantum field 
theories [14 , and projected entangled pair state (PEPS) 
for higher dimensional systems p!5HT7] . 

The MPS also appears to be the exact ground states of 
certain spin models. For example, the valence-bond solid 
(VBS) ground states of the Affleck-Kennedy-Lieb-Tasaki 
(AKLT) models [18 are matrix product states. They 
provide a clear physical picture to the Haldane gap phe- 
nomena [19 and shed light on their "nearby" integer-spin 
Heisenberg antiferromagnets [20 . In condensed matter 
physics, the Hamiltonians usually arise with two-body 
interactions and SU(2) symmetry since they are relevant 
to describe realistic materials. In order to study such 
systems, a method has been suggested in Ref. [21j to 
construct the SU(2) symmetric two-body parent Hamil- 



tonians for MPS. However, when starting with the Hamil- 
tonians, in principle it can be extremely hard to find their 
matrix product ground states [22 . 

The purpose of this paper is to investigate a real-space 
renormalization group and its applications in a system- 
atical search for exactly solvable quantum spin chains. 
The present approach complements the parent Hamilto- 
nian method in Ref. [21 , such that one can start from the 
Hamiltonians and search for exactly solvable ones. We 
first briefly review the basics of real-space renormaliza- 
tion and its extension to systems with SU(2) symmetry. 
The presence of symmetry allows us to design a simple 
exact renormalization scheme. By using this method, we 
study quantum spin chains with SU(2) symmetry and 
nearest-neighbor interactions. For 6* < 4, we provide 
complete solutions for the models which are frustration- 
free for two neighboring spins. Moreover, we also provide 
a new MPS solution of 5* = 5 which was not previously 
known. We discuss these exact solutions by dividing 
them into two different classes, whose ground states are 
matrix product states and partially ferromagnetic states, 
respectively. 



II. REAL-SPACE EXACT RENORMALIZATION 

Let us consider a chain with N local (i-dimensional 
Hilbert spaces H, that we can assume local spins. We 
denote by \M) G H diU orthonormal basis in H. And let 
us also consider a translationally invariant Hamiltonian 
H = ^. h- containing local interaction terms acting on 
contiguous k sites. We can assume positive semidefinite 
interaction /li > 0, since they can always be achieved by 
shifting local energy level. 

Let us now briefly explain the real-space renormal- 
ization process. We start by coupling the first two 
spins, whose Hilbert space H ^ H ^ I-L2 is mapped 
into a Hilbert space I-L2 which has in general a dimension 
D2 < (P. The criteria followed to perform this reduction 
is to conserve only the low energy states of the Hamilto- 
nian. In general, the method works by finding the map- 
pings ^[*] : Hi-i ^H ^ Hi which carry out this process. 
We continue this renormalization procedure until reach- 
ing the end of the chain and getting an orthonormal basis 
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Let us show the real-space renormahzation process 
from the {i — l)-th spin to the i-th spin, which can be 
written in a basis as |T0] 
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where the input state |Qf[i — 1]) G Hi_i, the output state 
\P[i]) G Hi, and the Kraus operator A^^^^'s are Di-i xDi 
matrices satisfying isometry condition ^^ ^[^]t^[^] = 
1. Here we define I^o = 1 so that the Kraus operator 
^[^i] for the first spin can be viewed as a row vector. 

Equation (IT]) shows the real-space renormahzation re- 
sults in an orthonormal basis |x) (x = 1 ' • ' Dn) with a 
matrix product form (See Fig. la) 

|x)= Yl {A^^'^A^^'^---A^^''\\MuM2---Mn) 

Ml... Mat 

(2) 
where D = max^ Di is called the bond dimension of the 
matrix product. In the DMRG algorithm, these matrix 
product states are used variationally to find the best ap- 
proximation of the low energy sector of the ID systems. 

In this work, we are interested in special models such 
that the states |x) exactly span the ground-state sub- 
space in the thermodynamic limit. The specifications 
about the thermodynamic limit comes from the fact that 
every state for N sites can be written by means of a ma- 
trix product ansatz given in Eq. ([2| by taking D > d^ . 
However, we seek for models for which an exact renor- 
mahzation can be performed for arbitrarily long chains. 
In other words, the ground states of these models can be 
solved rigorously through the real-space renormahzation, 
and the truncation induced by the Kraus operators does 
not harm. 

Practically, since /li > 0, this search can be achieved 
if the Kraus operators for each spin can be adjusted step 
by step in the renormahzation group to fulfill 

TT{p}hi) = 0, Vi = 1 • • • TV and Vx = 1 • • • l^AT (3) 

where pf = Trenv[|x)(x|] is the reduced density matrix 
for k spins. The above condition leads to iJ|x) = 0, which 
means that the vectors |x) are the ground states of H^ 
because H >0. Such Hamiltonians are called frustration- 
free Hamiltonians since their ground states minimize en- 
ergy locally. For instance, it is well-known that the fer- 
romagnetic Heisenberg chain is a typical frustration-free 
model in which all the spins tends to align in parallel 
to gain energy. Recently, the frustration-free Hamiltoni- 
ans have been reformulated as quantum /c-SAT problems 
and attract considerable interests in quantum informa- 
tion community |23ti26] . 

For generic models, this renormahzation procedure will 
terminate after blocking a number of spins due to intrin- 
sic frustrations. To find the exactly solvable model, the 
first possibility is that the renormahzation group reaches 
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FIG. 1: Real space renormalization. (a) The real-space 
renormalization group yields matrix product states, (b) In pe- 
riodic boundary condition, the translationally invariant MPS 
are constructed from the fixed point Kraus operators. 



a fixed point. Then, the ground state of the Hamilto- 
nian in periodic boundary condition can be written as a 
translationally invariant MPS (See Fig. lb) 

1^) = IZ Tr(A[^^U[^2] • • • A[^^])|Mi, M2 • • • Mat) 

Ml... Mat 

where the Kraus operators A^-^^ are the converged Dx D 
matrices at the fixed point. We discuss these fixed point 
MPS solutions in Sec. HI B. Another possibility is that, 
for some models, the number of states dimHk that we 
should keep, increases when gathering more spins. Even 
though there is no renormalization fixed point, we find 
that it is still possible to obtain the ground states exactly 
if dimHk increases in a controllable way. We illustrate 
this point in Sec. Ill C, when discussing the partially 
ferromagnetic states. 



III. QUANTUM SPIN CHAINS WITH SU(2) 
SYMMETRY 



In this section, we adapt the real-space exact renormal- 
ization to SU(2) symmetric quantum spin chains with 
nearest-neighbor interactions. Therefore, let us start 
by explaining some details about the SU(2) symmetric 



Hamiltonians. The most general SU(2) symmetric trans- 
lationally invariant spin-^S Hamiltonian with nearest- 
neighbor interactions can be expressed as 



25 



H = ^Y.an{S^ ■ S^^i^ ^ aol. 



(4) 



n=l 



The study of these SU(2) symmetric models has a long 
history in condensed matter physics. It was known that 
some of these models can be solved by Bethe Ansatz 
method and such models are fully classified by solutions 
of Yang- Baxter equations [27 . 

We want to identify the frustration-free models in Eq. 
Q and find their ground states through real-space exact 
renormalization. However, it is convenient to use projec- 
tors instead of spin operators, so we use the transforma- 
tion 

{s,-s,^ir = E [lsT{ST^i)-s{s^i)rPsAhi^i) 

St=0 

where P^^ is a projector onto total-spin St states of the 
two spins. By shifting the local energy levels, we can 
always rewrite the Hamiltonian Q as a sum of projectors 



H = ^Y1 JsTPsAhi^l) 



(5) 



St^K 



with coupling constants Js^ > and /C C [0, 25*] is a 
set specifying the choice of projector(s) as local interac- 
tions. Since the local interactions in Eq. ^ are a sum 
of projectors, we have H >0. 

Let us remark that, as the physical representation is 
irreducible and we restrict to nearest-neighbor interac- 
tions, the exact value of Jsj. is not important whenever 
the Hamiltonian is frustration- free. 

From the projector Hamiltonian (Isl), it is still not clear 
how to properly choose, if possible, the set /C to make 
the Hamiltonian frustration- free. However, as we restrict 
ourselves to frustration-free models with two neighboring 
spins, we can provide a complete list by taking advantage 
of the renormalization group. 



A. Exact renormalization Avith SU(2) symmetry 

In this subsection, we explain how to make use of the 
SU(2) symmetry in the exact renormalization scheme. 
This particularizes the real-space renormalization in Eq. 
^ such that both the input and output states form rep- 
resentations of the symmetry group, which ensures the 
symmetry is preserved in each renormalization step. The 
method shown here is a three-step process. 

Equation ^ can be promoted to an SU(2) adapted 

basis mm 



\j,hm,)= V VA^-i 
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\jatama)\SM) (6) 



where the j's denote the SU(2) representations (total- 
spin quantum number), the t's distinguish the degenerate 
states within the same j, and the tti's are the magnetic 
quantum numbers associated with j. 

The first step of the process consists of splitting the 
Kraus operators into two terms by means of Wigner- 
Eckart theorem (See Fig. 2a) as 



ASM] 
jatamajbtbrrib 



^ Tj^t^,j,u{Jama,SM\jkmb) (7) 



where the indices jata^jbh keep track of the representa- 
tions of the input and output states. The first term is 
a real matrix T denoting the weights of different input 
states in each output states. We call this matrix "weight 
matrix" . Let us remark that the weight matrix does not 
depend on the magnetic quantum numbers. The second 
term is the Clebsch-Gordan coefficient {ja'^ai SM\jijmij), 
corresponding to the representation fusion ja ^ S ^ j^. 
To ensure that the output states always form an or- 
thonormal basis, the weight matrix must fulfill 



'-jatajhth 



unless \ja - S\<jb<ja^ S (8) 



/ J -'-j at a, jhth-'- 3 gig jht'^ ~ ^th,t'^ 



(9) 



jgtg 



for every J5. The first constraint is related to SU(2) fu- 
sion rules. The second constraint means the columns of 
Tj^tajbtb corresponding to the same J5 but different t^, 
are orthonormal vectors, which guarantees the isometry 
condition ^^ At'^'^^'f'At'^'^] = 1 for the Kraus opera- 
tors. 

The advantage of this representation for the Kraus op- 
erators is that it allows us to design an elegant way to 
perform the exact renormalization group, which is the 
second step of the method. Let us consider two neighbor- 
ing spins (See Fig. 2b). The renormalization process con- 
sists of two sequential representation fusions ja^ S ^ jb 
and jb^ S ^ jc As a result, we obtain the orthonormal 
basis 



/ ^ -^jgtgjbtb 



Ijctcrric) = "^ Yl 

M1M2 jatgiria jbtbirib 

X {jarria, SMi \jbmb){jbmb, SM2 \jcmc) 
xT;^tb,j.tJjatama)\SMr)\SM2). (10) 

where the weight matrices T and T' for these two spins 
can be different in general. Alternately, the renormal- 
ization process Eq. (10) can be done by fusion of 



the two physical spins to their coupled representations 
S ^ S ^ St first and then ja ^ St ^ jc In the latter 
fusion sequence, we obtain the same basis 



Ijctcrric 



Yl Yl ^utg.lru{Jama.STMT\jcmc 

jatgrua StMt 



Jgtg ij ci^c 

Ajatama)\STMT) 



(11) 



where \StMt) is the coupled basis of two physical spins. 
The two different fusion channels are unitarily related 
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JataHl, 
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jaiiia Jb,ma+Mi jcin^ 



ZT^JaSJb 
-t^SJcSx • 

St 

(c) 



S,M2 



St,Mi+M2 



the orthogonal complement of /C. This simple test verifies 
the output states are the zero-energy local ground states 
of the projector Hamiltonian ([s]). 

The renormalization group follows naturally as T^^^ -^ 
T[^] -^ ■ - ■ under the renormalization condition Eq. (13) 



and the constraints Eq. ([8| and Eq. ([9|. Let us describe 
how to deal with these requirements simultaneously. Let 
us suppose that we already know the weight matrix T^*"^] 
and the goal is to calculate T^*]. After taking the square 
of Eq. (13) and summing over ja^ta and St ^ /C, we 
obtain 



Z^ Z^ ^j'A 
jbtb ^btb 



M 



i]j. 



T- 



,3 etc Jit'^^Jbtb jbtb J etc 



= 



(14) 



where the positive semidefinite real Hermitian matrix 
J\ylb]jc [^ given by 



M 



i]jc 

Jbt'jbtb 



STeKjata 



T^JaSj'^ j-ijaSjb I 
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tajiti SStJc SStJc jatajbtb' 



FIG. 2: Real space renormalization Avith SU{2). (a) An 

isometry A^^'^^ is decomposed as a matrix T and a fusion of 
angular momentums. (b) The renormalization of two spins 
can be done in two successive steps, (c) The same input and 
output states with different intermediate fusion channels for 
two spins can be related by the F-symbol. 



by the recoupling F-symbol (See Fig. 2c) defined by 

Fs€'jc = {Jb{jaS),S;jcmc\ja.ST{SS);jcmc). By using 
Wigner's 6-j symbol, this F-symbol can be expressed as 

EU 



F^tl = i-iy^^'Wi^ji^ + msr + 1) {^^ I j^ 



jja S jb \ 

\S jc St}' 



By substituting this in Eq. ( 10 ) and comparing with Eq. 
(11), we obtain 



TdSt _ \^ rp T^jaSjb rp/ 

^jatajctc ~ Z^-^Jatadbtb-^SSTJc-^JbtbJctc' 
jbtb 



(12) 



According to Eq. (11), the output states Ijc^c^^c) only 
keep the local ground states of the Hamiltonian ([5| if 



R 



Jata }Jctc 



= ^ja,ta and St & fC 



(13) 



This equation relates the weight matrices of two spins 
in Eq. (12) and plays an important role in our exact 



renormalization group method. 



The third step of our method is to use Eq. ( 13 ) to carry 



out the renormalization group process for the whole spin 
chain. Hereafter we use T^*] to denote the weight matrix 
at site i. Let us start the renormalization from the first 
two spins. By taking the first physical spin S as the input 
representation, we have the initial condition T^J = 1 with 



ji = S. According to Eq. (12) and Eq. (13), we obtain 

[21 — ~ 

ii 12 ~ ^' ^^^ output representation J2 G /C, where /C is 



For every possible jc from jb^ S^ we calculate the kernel 
of A^[*]-^^, which give us the weight matrix T^*]. If A^t*]-^^ 
does not have kernel vectors satisfying Eq. (Isl), the cor- 
responding output representation jc must be discarded. 
If the kernel of M.^^'^^'' has dimension larger than 1, the 
index tc is used to tag the orthonormal kernel vectors 
for such jc. Thus, the kernel vectors of A^t*]-^^ consti- 
tute the columns of T^*] and the column indices j^^c 
of TjJ^^j^^^ denote the output representations. One can 
straightforwardly show that the resulting weight matrix 
T[*] satisfies the renormalization condition Eq. (Il3| and 
the orthonormal constraint Eq. ([9|, because AT^^-^'^ is 
positive semidefinite and Hermitian. 



B. Matrix product states 

In this subsection, we discuss the models which have 
a renormalization fixed point and then, MPS as ground 
states. In our present exact renormalization scheme, the 
renormalization fixed point means that the output rep- 
resentations does not change when adding new spins and 
T[*] converges to a site-independent matrix. 

Let us start by introducing two relevant concepts about 
MPS - injectivity and symmetry. We begin with the 
definition of injectivity [37 

Definition 1 (Injectivity) Let \a) G C^ be an or- 
thonormal basis and {A^^^f^^-^ the D x D Kraus op- 
erators defining a translationally invariant MPS. And let 
us consider the D'^ states for L sites defined as 

l^i?)= E H^'^^l---^''^"l|/^)|Mi---Mi) (15) 

Mi---Ml 

Then, we say that the MPS is injective (see Fig. 3) if 
there exists a finite L such that the vector space spanned 
by the vectors in Eq. (15) has dimension D^ . In other 




FIG. 3: Injectivity. A MPS is injective with injectivity 
length Lo if for the MPS constructed for Lq sites, different 
boundary conditions (hnearly independent), represented by 
the cones in the figure, give rise to different states (hnearly 
independent), and this does not happen for Lo — 1 spins. 
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FIG. 4: Symmetry. The unitary u applied on the physical 
level is reflected in the virtual level as pair of unitaries [/, as 
shown in Ref. i2L 



words, different boundary conditions turn into different 
states. The injectivity length Lq is defined by the minimal 
number of sites for which injectivity is reached. 

The interest of this definition comes from Ref [5| 171 , 
where it is proven that injectivity is the necessary and 
sufficient condition for the existence of a parent Hamil- 
tonian which has the MPS as a unique ground state with 
a non-trivial spectral gap above. 

The other relevant result that we would like to recall 
here is the construction of translationally invariant MPS 
which are locally invariant under some symmetry group 
G. The following theorem provides the necessary and 
sufficient conditions [21] for that 

Theorem 2 (Symmetry) Let |^) G (C^)^^ be trans- 
lationally invariant MPS defined by the Kraus operators 
I^M}^^^^ and let u and U be two representations of a 
finite or a compact Lie group G. Then, \^) is invariant 



under G in the sense of u® \^) 
(see Fig. 4) 



iON 



^) if and only if 
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M' 



Um,M 



,A^M']^^i9jjj^[M]jj^ 



(16) 



Here we call that u and U are the physical and virtual 
spin representations, respectively. Once these representa- 
tions are fixed, the Kraus operators can be constructed by 



means of the Clebsch-Gordan coefficients together with 
a weight matrix. In our present SU(2) case, the Kraus 
operators are exactly given by the decomposition in Eq. 

0- 

Let us now remind the previous results about the MPS 
solutions for Hamiltonian ([5|. The best-known models 
belong to the AKLT family [18 , which are defined by 
/C = {6' + 1, 5' + 2, . . . , 2S'} for integer-spin S. The MPS 
of spin-S* AKLT model have a VBS picture with irre- 
ducible virtual spin-5'/2 representation. Another family 
of the models also have integer spin and the Hamilto- 
nians are defined by /C = {2, 4, . . . , 26'} [32l |32, which 
we call SO (26 + 1) symmetric family. For the 6 = 2 
model of this family, the MPS have irreducible virtual 
spin-3/2 representations [32 , which is equivalent to the 
SO (5) symmetric MPS in a two-leg electronic ladder [34]. 
For 6 > 3 cases, the properties of the corresponding MPS 
are less clear, even though their explicit wave functions 
were found [32 . 

Now we turn to our results obtained by the exact renor- 
malization group. For the Hamiltonian (Isl), we check all 
possible /C [35 and then provide a complete list of fixed 
point MPS solutions for 6 < 4, and a new solution for 
6 = 5. All these solutions are integer-spin models [36], 
which are summarized in Table [TJ For 6 < 4, we conclude 
that there is no solution other than the above two fami- 
lies. For 6 = 5, we find a new model, whose Hamiltonian 
is given by /C = {3, 7, 8, 9, 10} and the ground state has 
a VBS picture with irreducible virtual spin-3 representa- 
tions. 

For the SO(26 + 1) family with 6 > 3, the exact 
renormalization group provides us a more comprehensi- 
ble physical picture, which can be viewed as generalized 
VBS with reducible virtual spin representations. In Ta- 
ble |l| we also listed the minimal number of blocked spins 
to reach the fixed point representations. Since all these 
MPS are injective, this length scale is actually the injec- 
tivity length [37 . 

Let us explain these results with an explicit example 
in the SO (26 + 1) family: the spin-3 model with /C = 
{2,4,6}. Through the exact renormalization group, we 
can observe that the output states reach the fixed point 
representation 0©0©1©2©3©3©3©4©5©6 
after blocking 6 spins. To obtain the MPS, we do not 
really need to calculate the fixed point Kraus operators 
by the renormalization group. According to Theorem 2, 
the fixed point representations allow us to construct this 
MPS directly. For the present example, the fixed point 
representations give an important hint that the MPS has 
a VBS picture (See Fig. 5a) with SU(2) reducible virtual 
spin representation © 3, which is quite different from 
the traditional VBS states with irreducible virtual spin 
representations, like AKLT states [18 or their extensions 
[5.. 

With a chain beyond the injectivity length Lq = 6, 
the tensor product of two © 3 representations at the 
two boundaries yields the observed fixed point represen- 
tation in the renormalization group. For open boundary 
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FIG. 5: Ground-state physical picture, (a) The fixed 
point type MPS solutions have a VBS picture. Each dot de- 
notes a virtual spin representation. The wavy lines represent 
the valence-bond singlets between virtual spins and the circles 
indicate the projection of two virtual spins onto physical spin 
representations. In the exact renormalization calculations, 
the fixed point representations are from the tensor product of 
two virtual spins (edge states), (b) The partially ferromag- 
netic states have a magnetization plateau. The arrows denote 
a fully polarized virtual spin-(5' — 1) in a spin-S* partially fer- 
romagnetic state. 



Spin 


Set /C Virtual spin Lq 


1 


{2} 1/2 2 


2 


{3,4} 1 2 


2 


{2,4} 3/2 4 


3 


{4,5,6} 3/2 2 


3 


{2,4,6} 00 3 6 


4 


{5,6,7,8} 2 2 


4 


{2,4,6,8} 20 5 8 


5 


{6,7,8,9,10} 5/2 2 


5 


{2,4,6,8,10} 5/2 9/2 15/2 10 


5 


{3,7,8,9,10} 3 4 



TABLE I: Models with 5'[/(2)-invariance, nearest-neighbour 
interactions and matrix product ground states. Lq is the in- 
jectivity length. 



conditions, in thermodynamic limit, the unpaired rep- 
resentations © 3 at the two edges are asymptotically 
free and become well-defined edge states. For periodic 
boundary conditions, all virtual spin representations are 
contracted into SU(2) singlets with neighboring sites and 
therefore the MPS is a global spin singlet. 

The renormalization group analysis has also been car- 
ried out for other models in 80(2/5 + 1) family From 
Table [l| one can see that, for S* > 3, their matrix prod- 
uct ground states have reducible virtual spin representa- 
tions, which directly correspond to the edge states in an 
open chain. This provides more complete understanding 
of these systems. For all MPS in Tab. |Tj we present their 
explicit Kraus operators in Appendix A. 

Let us make a remark about these exactly solvable 
models. All their fixed point MPS ground states have 
exponentially decaying correlations and there is an en- 
ergy gap above the ground states, since they are injec- 
tive. However, the different virtual spin representations 
(edge states) show that these MPSs belong to different 



quantum phases of matter. Therefore, once a new Hamil- 
tonian H = {1 — x)Hi -\- xH2 is constructed from two 
solvable models Hi and H2 in Table [T| with the same spin 
S^ at least one quantum phase transition is expected to 
occur when tuning x from to 1. Since both MPS ground 
states for Hi and H2 preserve SU(2) symmetry, the lo- 
cal order parameter description breaks down and uncon- 
ventional quantum phase transitions may emerge. Very 
recently, this idea has been exploited to study the pos- 
sibility of a topological quantum phase transition in an 
5' = 2 chain |3H]. 



C. Partially ferromagnetic states 

In this subsection, we discuss another class of models 
which do not have renormalization fixed points but still 
can be solved exactly. The ground states of these models 
are partially ferromagnetic states. 

This family includes both semi-integer spin models 
and integer-spin models. The Hamiltonian is defined by 
/C = {0, 1, ... 25' - 4, 2S} and the physical spin S > 5/2. 
Their ground states are partially ferromagnetic states 
with a magnetization plateau (Sf) = S — 1. We also 
have found a physical picture (See Fig. 5b) for these 
states with partial magnetization: We prepare a spin-1 
AKLT-type VBS state with virtual spin-1/2 and a spin- 
{S — 1) maximally polarized ferromagnetic state. In each 
site, we recover the physical spin-^ Hilbert space by 
(5 — 1) (g) 1 -^ 5, which is achieved by applying local 
projections. 

Let us consider a typical example - the spin-5/2 model 
with /C = {0,1,5}. For a block of Nq spins, the 
AKLT part contributes representations 1 and the 
polarized ferromagnetic part contributes representation 
3Nq/2. Thus, the total spin of the A/'o-spin block is given 
by the tensor product of representations from these two 
parts 



(O©1)0 



37Vo ^ 37Vo 



-1)( 



3A^o 

) ( 

2 



3A^o 

) ( 

2 



.(?f + 1) (17) 



For two spins {Nq = 2), the allowed representations are 
2 © 3 © 3 © 4 and can not reach /C = {0, 1, 5}, which 
means that the partially ferromagnetic state is the zero- 
energy ground state of the projector Hamiltonian. In the 
exact renormalization process, we found that the four 
output representations in Eq. (17) are the only output 
representations for Nq > 6. By adding one additional 
spin, the total spin of the four representations is increased 
by 3/2. These observations actually strongly suggest the 
partially ferromagnetic picture of the ground state. 

One may ask why this class starts with S = 5/2 rather 
than S* = 2. The reason is the following: For spin-2 
model JC = {0,4}, the renormalization group shows that 
the number of output representations does not saturate, 
which means that the partially ferromagnetic state is not 
the only ground state of the Hamiltonian. 



Compared to the fixed point MPS solutions in Sec. Ill 
B, the partially ferromagnetic states have a long range 
order and thus break the SU(2) symmetry. According to 
Goldstone theorem, we expect gapless spin wave excita- 
tions above the ground state, which is quite different from 
the gapped fixed point MPS with exponentially decaying 
correlations. 



IV. CONCLUSIONS AND PERSPECTIVES 

We have introduced a real-space exact renormaliza- 
tion group adapted to the SU(2) symmetry, which is well 
suited for finding exactly solvable quantum spin Hamil- 
tonians with nearest-neighbor interactions. 

The list of solutions can be divided into two classes 
according to the renormalization group behavior. In the 
first class, the models are quantum integer- spin chains 
with renormalization fixed points and matrix product 
ground states. For S* < 4, we show that the AKLT 
family and the S0(2S' + 1) family exhaust all possible 
solutions. In the SO(26' + 1) family, the renormalization 
group provides a natural explanation for the edge states 
of the MPS by providing a generalized VBS picture with 
reducible virtual spin representation. Furthermore, we 
obtain a new solvable model for 6* = 5 beyond the exist- 
ing families. In the second class, the models have par- 
tially ferromagnetic ground states with a magnetization 
plateau. This solvable family exists for S > 5/2 and 
contains both integer spin and semi-integer spin models. 
The partially ferromagnetic ground states have gapless 
spin- wave excitations, which are quite different from the 
gapped MPS in the first class. 

Beyond the present work, it would be quite interesting 
to generalize the method to spin chains beyond nearest- 
neighbor interactions and models in higher dimensions, 
especially adapted to PEPS formalism. Furthermore, the 
method may be used to explain an open question by 
Ostlund and Rommer [10] about which representations 
must be introduced in the renormalization group with 
SU(2) symmetry. 

Finally, we also expect a natural extension of the 
present exact renormalization formalism to quantum spin 
chains with other symmetry groups [39 . 
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Appendix A: Kraus operators of the matrix product 

states 

In this Appendix, we explicitly present the Kraus oper- 
ators needed for the definition of the MPS in Table HI As 
we mentioned, the Kraus operators with SU(2) symme- 
try are parametrized by Eq. (I?]), which requires both the 
set V containing the SU(2) virtual spin representations 
and the weight matrix T. 

For irreducible virtual spin representations, the set V 
contains a single representation ja and therefore T = 1. 
In this case, the Kraus operators are simply the Clebsch- 
Gordan coefficients 



a[S,M] 
jamajarrib 



{jamais MlJaTTlb) 



(Al) 



For reducible virtual spin representations, the set V has 
multiple SU(2) representations and the weight matrix T 
is necessary. The Kraus operators are given by 



i[SM] 



Jama,Jbmb 



^{jama,SM\jbmb). (A2) 



where the index t is suppressed because no degeneracy 
occurs in V for our models. We use a convention to de- 
fine the matrix T such that the row and the column in- 
dices jadb are arranged in an incremental order. For 
instance, the S* = 3 model with /C = {2, 4, 6} has virtual 
representation 3 and 



^0,0 ^0,3 
X3,0 X3,3 



For the 5 = 4 model with /C 
representation 2 5 and 





(A3) 



{2, 4, 6, 8}, we have virtual 



(A4) 



For the 6' = 5 model with /C = {2,4,6,8,10}, we have 
virtual representation 5/2 9/2 15/2 and 



(A5) 



It is straightforward to show that these solutions satisfy 
Eqs. ([8|, ([9|, and ([l3). 
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